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Abstract 

We have studied strange quark stars in the framework of the MIT bag model, allowing the bag 
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parameter B to depend on the density of the medium. We have also studied the effect of Cooper 
pairing among quarks, on the stellar structure. Comparison of these two effects shows that the 
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former is generally more significant. We studied the resulting equation of state of the quark matter, 



£Sj . stellar mass-radius relation, mass-central-density relation, radius-central-density relation, and the 

> 

variation of the density as a function of the distance from the centre of the star. We found that 

t>' 

the density-dependent B allows stars with larger masses and radii, due to stiffening of the equation 

G\ ; 

of state. Interestingly, certain stellar configurations are found to be possible only if B depends on 



the density. We have also studied the effect of variation of the superconducting gap parameter on 



Oh' our results. 
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I. INTRODUCTION 



The conjecture that the true ground state of quantum chromodynamics (QCD) may be 
the strange quark matter (SQM) was first explicitly stated by Witten flj, though various 
forms of quark matter had been considered earlier 0, [J . The key idea is that the energy per 
baryon of SQM could be less than that of even 56 Fe, the most stable nucleus, making SQM 
comparatively more stable. Farhi and Jaffe jj| studied SQM in the framework of the MIT 
bag model |5[, for various values of the strange quark mass (m s ) and the bag constant (B). 
They found a "window of stability" in the m s -B plane inside which SQM was stable. The 
astrophysical implication of the stability of SQM, namely the possible existence of strange 

n n 

(quark) stars, was studied by Haensel et al. p and Alcock et al. [7]. They found that 
very compact strange stars are possible, which are much smaller than the normal neutron 
stars, the reason being that the strange stars are self-bound and generically different from 
the gravitationally bound stars. Of course, the actual existence of strange stars can only be 
confirmed by observational studies of pulsars and pulsar-like compact stars, including X-ray 
pulsars, X-ray bursters, anomalous X-ray pulsars, soft 7-ray repeaters and isolated neutron 
stars. 

The SQM hypothesis has been around for almost two decades. Recent theoretical progress 
in our understanding of the cold and dense QCD matter has renewed the interest in this 
hypothesis: At low temperatures and sufficiently high baryon number densities, it is thought 
that quarks near the Fermi surface form Cooper pairs and new condensates develop. ^This 
is expected to lead to a colour-superconducting phase in the QCD phase diagram 



Among the various colour-superconducting phases, the colour-flavour-locked (CFL) phase 



10 



11 



12j is the one in which quarks of all three colours and all three flavours (it, d, s) 
pair. Electrical neutrality of bulk CFL matter is ensured due to the fact that this phase 
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favours equal number densities of u, d and s quarks, despite their unequal masses 
This also means that electrons are absent in the bulk CFL matter. Strange stars made of 
CFL matter seem to be favoured over those made of another colour-superconducting phase, 
namely the 2SC phase jig . These theoretical developments have important implications for 
the existence of quark stars and neutron stars with quark cores. 

Other reasons for the revival of interest in this area are the recent claims of observation 
of both mass and radius of compact stars jlfil . 0|. Although these claims are somewhat 
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controversial Q|, they have further stimulated theoretical activity in this area. 

The stability of SQM in the light of the above theoretical developments, and the structure 
of strange stars have recently been investigated by Lugones and Horvath 

Q,Q- While 

constructing the equation of state (EOS) of the CFL matter, they introduced the vacuum 
energy by means of the phenomenological bag constant B. However, in view of the high 
densities encountered in neutron stars and quark stars, the question whether the bag constant 
should be held fixed at its free-space value or whether its variation with the density of the 
medium should be taken into account becomes very relevant. In this work, we investigate 
the effect of a density-dependent B, on the EOS and the resulting structure of the CFL 
strange stars. 

The paper is organized as follows. In Sec. II we discuss various models of B(p) where p 
is the baryon number density of the medium. In Sec. Ill we consider the CFL matter and 
its EOS. Our numerical results for the resulting stellar structure are given in Sec. IV. In 
Sec. V we present our conclusions. 

II. DENSITY-DEPENDENT BAG PARAMETER 

Of the many models of the nucleon that have been constructed so far, one of the most 
useful is the MIT Bag Model The nucleon is assumed to be a "bubble" of perturbative 
vacuum immersed in the non-perturbative or true vacuum, and quarks are confined to the 
bubble by means of a net inward pressure B exerted on it by the surrounding vacuum. The 
bag model has also been used to study the thermodynamics of the deconfinement phase 
transition. As the temperature rises above the deconfinement temperature, there is no 
difference between the two vacua, and the net inward pressure B must vanish. In other 
words, B must be viewed as a temperature-dependent quantity [21]. The analogy between 
the bag constant B and the condensation ener gy i n the Nambu-Jona-Lasinio (NJL) model 



at finite temperature was pointed out in Ref. [22j. It was shown that as the temperature 
rises, this condensation energy decreases and goes to zero at the transition temperature. 
This provides a model of the temperature-dependent B. Density dependence of B can be 
motivated in an analogous manner: Since the deconfinement and the concomitant vanishing 
of B can also be brought about b y ra ising the baryon number density p, B should be treated 
clS db density-dependent quantity 
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There have been a number of attempts in the literature to evaluate the density dependence 
of B. The results, however, are model dependent, and no consensus seems to have emerged. 
In this work, we made use of the results of the following three studies published recently: 
(a) Liu et al. j^J extended the global colour symmetry model [2^ to finite quark number 
chemical potentials /1, in order to derive the variation of B with p and p. (b) Burgio et al. 
j3| used the CERN SPS data on heavy-ion collisions to justify and determine the density 



dependence of B. (c) Aguirre 27|, |28j used the NJL model to study the modification of the 
QCD vacuum with increasing baryonic density and to extract the medium dependence of 
B. 

(a) Since Liu et al. have not presented their B(p) or B(p) in a parametric form, we have 
fitted their results with a suitable analytic expression with two parameters a\ and a^. 

B(p)/B(0) = exp[-(a!X 2 + a 2 x)\, (1) 

where x = p/po is the normalized baryon number density, po is the baryon number density 
of the ordinary nuclear matter, a\ = 0.0125657, and a 2 = 0.29522. They took -B(O) = 114 
MeV/fm 3 =(172 MeV) 4 . (b) Burgio et al. have presented their B(p) in a parametric form: 

B{p) = B as + (B - B as ) exp(-(3x 2 ), (2) 

where B as = 38 MeV/fm 3 , B = 200 MeV/fm 3 = (198 MeV) 4 , (3 = 0.14, and x is as defined 
above, (c) Aguirre [28] has calculated B(p) for a symmetric uds quark matter relevant for 
a CFL quark star. His B(p), in MeV/fm 3 , is given by 

B(p) = a + hx + b 2 x 2 + b 3 x 3 + b A x 4 + b 5 x 5 , x < 9 

= exp(-a(x - 9)), x>9 (3) 

where x is as defined earlier. Further, a = 291.59096, 61 = -142.25581, b 2 = 39.29997, 63 = 
-6.04592, 64 = 0.46817, b 5 = -0.01421, a = 0.253470705, and (3 = 19.68764. Thus in this 
case B(0) = 291.59096 MeV/fm 3 = (217.6 MeV) 4 . 

We display, in Fig. 1, B(p) vs p corresponding to Eqs. (l)-(3). These three models 
provide B(p) over a sufficiently wide range of densities expected to be seen inside quark 
stars. Secondly, the above three values of B(0) span more or less the whole range of values 
of the bag constant found in the literature. Finally, as is clear from Eqs. (l)-(3), they yield 
qualitatively different shapes for B(p): (1) and (3) have negative slopes at x — 0, while (2) 



has a vanishing slope; (1) and (3) vanish as x — > oo, while (2) tends to a positive constant. It 
will be interesting to investigate what effects these different shapes and magnitudes of B(p) 
have on the EOS and the stellar structure. We shall use these three models in subsequent 
sections; all our numerical results are labeled accordingly as (a), (b), or (c). 



III. THERMODYNAMICS OF THE CFL PHASE 



ree energy density Vt CFL of the CFL phase, at temperature T = 0, is obtained as 
13^ . We first consider the free energy density fij ree of an unpaired quark matter 
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in which all quarks that are "going to pair" have a common Fermi momentum v: the 
CFL pairing is most effective if the Fermi momenta of the u, d, and s flavours are the 
same. Although the nonzero strange quark mass tends to favour unequal number densities 
for the three flavours, the CFL pairing forces the three flavours to have the same Fermi 
momentum and hence the same number density, so long as m s is not too large. If A is 
the colour-superconductivity gap parameter and p is the quark number chemical potential, 
then the binding energy of the qq Cooper pairs can be included (to leading order in A/p) 
by subtracting 3A 2 p 2 /7i 2 from Qf ree . The vacuum energy density in the presence of a finite 
chemical potential p, is included by adding B(p). This gives 

3 

= ~2 f\p - »)P 2 dP + - 2 / V + ™ 2 s) 1/2 ~ l¥*P - 4 A V + BM, (4) 

sec which, however, treat B as a constant. The baryon number density p = 

—dQcFi/i^dfi) and particle number densities for individual flavours are given by 

+ 2AV) 1 

p = n u = n d = n s = -B (p), (5) 

where B'{p) = dB/dp. The common Fermi momentum v is chosen to minimize the free 
energy density with respect to a variation in v. This gives 

i 

3 

Energy density of the CFL matter is given by 

e = PiUi + 9-cfl = 3pp - P, (7) 



v = 2p-(p* + ^) . (6) 



since /i = (fi u + fid + Hs)/3 and the pressure P is given by P = —Qcfl- 

Now we obtain the equation of state of the CFL matter, in the standard form. Integrations 
in Eq. (4) can be performed in a straightforward manner and we get 
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; m 2 + ^2)3/2 _ ^ (m 2 + ^1/2 _ K ^-1 

4 8 8 m s 
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We are interested in the region where m s is small com par ed to the chemical potential fi. So 
it is sufficient to keep terms up to order m\ in Vtf ree To that end, first note that the 
common Fermi momentum v (Eq. (6)) is given by 
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where we have kept terms up to order m\ . This yields 

-3/i 4 3m 2 /i 2 l-121og(m s /2/i) 4 3 22 
" CFL = 7^ + "1^ 32^ ^ " ^ A ^ + (10) 

Substitution of v given by Eq. (9), in Eqs. (5) and (7) yields 

u 3 m 2 u m\ 2 a2 1 . 

7T 2 27T 2 87r 2 yU 7T 2 3 

9u 4 3m 2 u 2 11m 4 31og(m s /2u) 4 3 a2 2 . _,, . . . 

The pressure is given by 

3/i 4 3m^/i 2 1 - 12 log(m 8 /2/i) 4 3 2 2 
P= 4^-^ + 32^ ™* + ^ ~ B M- ( 13 ) 

Equations ()12|) and (fT3*j) give e and P in terms of /i as a parameter. The EOS in the 
standard form can be obtained from these two equations and it reads 

,t,/ s 6A 2 u 2 3m 2 u 2 1 + 61og(m s /2u) , . ,„ 

, = 3P + 4B(ai) - + -^f + 4 ; l 2 *V - ^V). (14) 

If only the first two terms on the rhs of Eq. (fH^ are kept then it reduces to the ideal gas 
EOS in the bag model. The occurrence of the pairing term stiffens the EOS because we get 
a higher pressure at the same energy density. A nonzero value of m s , on the other hand, 
tends to reduce the pressure leading to a softer EOS. 

If Eq. (14) is to be used to study the stellar structure, one needs to know B(fi) and its 
derivative B'(fi). In Sec. II, we have described three models of B(p). Consider first the 



model of Liu et al. j^J, denoted by (a) in Sec. II. Their formalism in fact first yields B(p), 
not B(p). They, however, do not present their numerical results for B{p). Instead, they use 
the relation p = 2/x 3 /(37r 2 ) to obtain B(p) and present numerical results for B(p). In order 
to extract their B(fi), we have to use exactly the same relation between p and /x. This is 
how we obtained B(p) in the case of the model (a). Next, consider the models (b) and (c) 
of Sec. II. To obtain B{p) from B(p) in these cases, we proceed as follows. Consider Eq. 
(11). It can be rewritten as 

/<">-!£• as) 

where f{p) denotes the sum of the first four terms on the rhs of Eq. (11). In the case of 
model (b), Eq. (2) can be inverted to get p(B) which can be substituted in Eq. (15). This 
gives a first-order differential equation between B and /x, which can be solved numerically. 
In the case of model (c), it is easier to eliminate B between Eqs. (3) and (15), and the 
resulting differential equation looks like 

(IT 

— = 3[f(ti)-p x]/(dB/dx). (16) 

In both the cases, the initial condition necessary to solve the differential equation is obtained 
from the fact that for sufficiently large pi, dBjdp, is nearly zero. We have checked the solutions 
of the above differential equations by substituting them in the original equations (Eqs. (15) 
and (2) in the case of model (b), and Eqs. (15) and (3) in the case of model (c)). 

Before we display our results for the EOS, Eq. (14), it is useful to see how the quantities 
B, B', e in Eq. (12), and P in Eq. (13) vary with /x. This is shown in Fig. 2. Consider first 
B(p). Its broad features are in agreement with those of B(p) seen in Fig. 1. Also note that 
in spite of the large differences among the values of B(p = 0) corresponding to the three 
models, the curves for B{p) in the region relevant to quark stars (/x ~ a few hundred MeV) 
do not differ drastically from each other. Consider next B'{p). It is either negative or zero 
depending on the value of /x. We have plotted \p,B'(p,)\ since it occurs in Eq. (12). This 
quantity obviously tends to vanish at very small and very large values of /x, and it displays 
a gentle peak at an intermediate value of /x. e{p) contains |/x-B'(/x)|, whereas P(p) does 
not. However, because of the relatively small magnitude of \p,B\p)\, the curves for e{p) 
and P{p) have similar shapes. Interestingly, the magnitudes of e(p) and P{p) are seen to 
be insensitive to the details of the three models of B(p). Finally, note that P{p) is positive 
only when /x exceeds a certain threshold value. 



Figure 3 displays the EOS of the CFL matter and unpaired quark matter, in the various 
models considered above. It is easy to understand the shapes of these curves by looking at 
P{p) and e(/i) in Fig. 2, and noting that Fig. 3 is obtained essentially by eliminating p. 
Some general trends emerge: (1) If B is treated as a constant independent of the density of 
the medium then the EOS appears to be linear; see the dashed lines in Fig. 3. On the other 
hand, if B depends on p, the EOS is slightly nonlinear. (2) The pairing term is seen to 
make the EOS stiffer, so does the density dependent B. However, the latter effect is seen to 
be more significant. This will be corroborated further when we present the stellar structure 
curves in the next section. 

We now discuss the issue of the "window of stability". This is the region in the m s -B 
plane, in which the SQM has energy per baryon (e/p) less than 939 MeV, and which was 
first discussed in Ref. |4j. Since both m s and B tend to increase e/p, this region obviously 
lies near small values of m s and B, in the first quadrant of the m s -B plane. It is bounded 
on the left by the vertical line B = 57 MeV/fm 3 =(145 MeV) 4 , and on the other side, by the 
e/p = 939 MeV contour. The boundary on the left arises because if B were smaller than 57 
MeV/fm 3 , i.e., if the inward pressure holding the quarks together in the nucleon were too 
small (see Sec. II), the nucleon and hence the nuclei would have dissolved into ud quark 
matter. Since that does not happen, B must not be smaller than 57 MeV/fm 3 . 

Recently, it was shown that the window of stability is enlarged considerably if the SQM 
undergoes a phase transition to a CFL quark matter jl^. The enlargement occurs because 
the e/p = 939 MeV contour gets shifted to the right. Let us understand this. To probe the 
absolute stability of the SQM, we need to set the external pressure P equal to zero. Then 
Eqs. (7) and (10) yield, respectively, 

e/p = 3/i, (17) 

and 

^ 1 + ^4aV + %) = 0, (18) 



47T 2 47T 2 7T 2 

where we have ignored the (9(m 4 ) term in Eq. (10), for simplicity. Ignoring the //-dependence 
of B as in ^| and eliminating p in the last two equations one gets e/p as a function of m s , 
B, and A. It is clear that a nonzero value of A shifts the e/p = 939 MeV contour to the 
right thereby enlarging the window of stability. This is essentially because of the opposite 
signs of the B- and A-terms in Eq. (18): The A-term tends to nullify the effect of a large B. 

8 



The larger the A, the larger the value of B one may have, while still maintaining e/p < 939 
MeV. This is clearly borne out by the results of 19]. 

Let us now discuss how a /i-dependent B would affect the window of stability. Earlier, 
when A was introduced, the functional dependence of e/p on m s and B was changed (see 
Eqs. (17)-(18)), and hence the contour e/p = 939 MeV got shifted. Now, whether B 
is constant or .-dependent, the new functional dependence remains unchanged. In other 
words, the e/p = 939 MeV contour in the present work is the same as in [19j . We shall 
revisit the issue of the window of stability in the next section, when we discuss Fig. 4. 

IV. STRUCTURE OF STRANGE STARS 

The equations of stellar structure, namely the Tolman-Oppenheimer-Volkoff (TOV) equa- 



tions 32, 33], give rise to a one-parameter family of stars corresponding to a particular 
equation of state: By specifying the central density as the parameter, one can numerically 
integrate the TOV equations, starting at the centre of the star and going radially outward 
until the pressure becomes zero which indicates that the surface of the star has been reached. 
This determines the radius R and the mass M of the star, for the specified central density. 
By choosing successively larger values of the central density, one can generate a sequence 
of stars of increasing mass. The sequence ends when any further increase in the central 
density leads to a star with a lower mass, which indicates an unstable configuration. This 
end-of-sequence star is the maximum-mass star for the sequence. For a given EOS, the mass- 
radius relationship is the most important property of the family of stars. This is because it 
can be compared directly with observational data in order to test and calibrate the theory. 
While both the mass and the radius are not known so far for any pulsar, a number of mass 
measurements exist and the data on the radius are getting better. 

We now present our numerical results on the structure of the strange stars. We set the 
strange quark mass m s = 150 MeV and the gap parameter A = or 100 MeV (except in 
Figs. 8 and 9 where we vary A between and 150 MeV). We recall from Ref. [19] that the 
stability of SQM requires that if m s = 150 MeV and A = 100 MeV, then the bag constant 
B should be less than about 117 MeV/fm 3 , and if m s = 150 MeV and A = 0, then B should 
be less than about 77 MeV/fm 3 . Now consider the three models of the density-dependent 
B, described in Sec. II. Since B is a decreasing function of the baryon number density p, 
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the above upper limits on the values of B necessarily imply lower limits on the values of p. 
In the numerical work, one has to ensure that these limits are satisfied, so that the resulting 
star is always within the window of stability. 

We now discuss Figs. 4-9. In all these figures, dashed curves correspond to the case when 
B is taken as a constant independent of the density of the medium, while the solid curves 
correspond to the case when B varies with p, and the panels (a), (b), (c) refer to the three 
models of B(p), described in Sec. II. 

Figure 4 shows the mass-radius sequences. Consider first Fig. 4(a). Model (a) has 
B(0) = 114 MeV/fm 3 which is well within the window of stability if A = 100 MeV (see 
above). This allows us to draw the two curves labelled 100. They show that a density- 
dependent B gives rise to larger and more massive stars. If A = 0, the above value of B(0) 
is outside the window. Hence there is no dashed curve labelled 0. However, the solid curve 
labelled is allowed: Recall that the density cannot be arbitrarily small since the pressure 
needs to be positive (Fig. 2). Indeed, for any point on this curve, the density throughout 
the star is found to be high enough that B(p) is everywhere less than 77 MeV/fm 3 (see 
above). Interestingly, this sequence of stars has become possible only because we considered 
a density- dependent B. Now consider Fig. 4(b). Model (b) has B(0) = 200 MeV/fm 3 
which is outside the window of stability for both A = and 100 MeV. Hence there are no 
dashed curves here. However, the solid curves labelled and 100, are allowed because for 
any point on these curves, the density throughout the star is found to be high enough that 
B remains within the window of stability. Note again that these sequences of stars have 
become possible only because we considered a density- dependent B. Similar remarks can 
be made about Fig. 4(c). Results in Fig. 4 seem to be robust — not very sensitive to the 
details of the three models of B(p) considered in this paper. 

The dot-dashed lines crossing the M-R curves in Fig. 4 correspond to the M/M© ~ 0.15 R 
(in Km) constraint obtained by Cottam et al. 3^ by using measurements of the red shift 
of spectral lines in X-ray bursts from EXO0748-676. In all Figs. 4(a-c) the crossing occurs 
at M/M e ~ 1.3-1.4. Thus the present calculation has been able to construct strange stars 
which are compatible with this constraint. 

It is a built-in feature of the TOV equations that the sequence of stars ends, whatever may 
be the EOS. The typical shape of the sequence shown in Fig. 4 (M oc R 3 approximately) is 
generic to compact stars that are self-bound. This is borne out by the fact that although 
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the EOS used here is quite different from that used in |3, |7| or |20|, the shapes of the curves 
in Fig. 4 are similar to those obtained in the earlier studies. 

Density-dependent B makes the EOS stiffer; see Fig. 3. A stiffer EOS can resist the 
collapse of a star to a black hole, to a greater extent. Hence stars with larger radii and 
masses result. What is true for an individual star holds also for a sequence of stars, which 
as a whole shifts toward larger mass and radius. This is evident in Fig. 4(a). 

Further, since the EOS is not strictly linear, simple scaling where the mass and radius 
scale as 1/ \f~B no longer holds jjy] . Actually when /i-dependent B is considered, the meaning 
of scaling as such is lost because as one moves from the surface of a star toward its centre, 
B keeps decreasing. That is, no fixed value can be assigned to it. 

Finally, in Fig. 4(a), note that the difference between the two curves labelled 100 is more 
significant than the difference between the two solid curves. Similar observations can be 
made about Figs. 5-7. This underscores the importance of the modification of B in a dense 
medium vis-a-vis the phenomenon of Cooper pairing among quarks. 

We plot mass M vs central energy density e c , in Fig. 5. As e c rises, so does the mass, 
until a plateau is reached which leads to the end of the sequence. The slope dM/de c is 
positive signifying stable configurations. Since the solid curves correspond to a stiffer EOS 
(Fig. 3), for a given value of the central energy density they represent more massive stars. 
Figure 6 shows the corresponding variation of the radius R with e c . 

In Fig. 7 we show, for the maximum-mass star, the variation of the density with the 
distance r from the centre of the star. In the case of the density- dependent B, the EOS is 
stiffer (Fig. 3) and the central density is lower. In other words, the solid curve labelled 100 
in Fig. 7(a) begins at a lower value of p/ p than the dashed curve. In spite of this, the star 
in the former case is more massive because it has a larger radius, and the large-r regions 
contribute to the mass in a dominant fashion. 

We display the variation of the maximum mass with the gap parameter A, in Fig. 8. The 
increase in mass with A is expected because as A increases the EOS becomes stiffer. Recall 
that in model (a), -B(O) = 114 MeV/fm 3 . For this point to be inside the window of stability, 
A has to be at least ~ 100 MeV. Hence in Fig. 8(a) the dashed curve starts near A = 100 
MeV. In Fig. 8(b) there is no dashed curve because in model (b), B(0) = 200 MeV/fm 3 
which is not inside the window if A < 150 MeV. A similar remark can be made about Fig. 
8(c). In Fig. 9 we plot the maximum radius vs A. Recent data hint at the possibility of 
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very compact stars ly, 17]. A smaller and denser object in comparison to the neutron stars 
will be a strong candidate for a strange star. 



V. CONCLUSIONS 



The fact that quarks normally exist in a confined state, does not automatically rule out 
the strange quark matter hypothesis. This hypothesis does not violate any physical principle 
and is also not ruled out by any observed phenomenon. The real difficulty in verifying the 
hypothesis lies in the fact that the underlying theory, namely the QCD at quark chemical 
potential (fi) of the order of several 100 MeV, is difficult to handle and does not yet allow 
a rigorous calculation. Unfortunately, the phenomenological models which one resorts to, 
give the energy per baryon of SQM very close to that of 56 Fe. In view of the uncertainties 
and the approximations involved, it is difficult to make a definitive statement. 

One way to make progress in this situation is to make the models more realistic. Taking 
the density dependence of B into account is a step in that direction. 

1. We have shown that this has a significant effect on the equation of state of SQM and 
hence on the resulting stellar structure curves. Specifically, the density dependence of 
B leads to a stiffer EOS which gives rise to stars with larger masses and radii. 

2. Our results for the M-R curves do not depend sensitively on the details of the three 
models of the density dependence of B considered in this paper. 

3. Effect of the newly discovered phenomenon of colour superconductivity on the struc- 
ture of compact stars has been studied by several authors. An important point which 
emerges from the present work is that the results depend more sensitively on the den- 
sity dependence of the bag parameter than on the phenomenon of Cooper pairing of 
the quarks. 

4. Another observation that we made is with regard to the 'window of stability' of SQM. 
We showed that the density dependence of B does not affect the window. 

5. We have been able to construct stars that are compatible with the constraint given in 
js^ . Moreover, their masses lie in the astrophysically plausible range. 
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6. Most importantly, certain stellar sequences which are ruled out by the requirement of 
the stability of SQM, if B is treated as a constant, are allowed if the density-dependence 
of B is taken into consideration. 

There are a number of ways in which this work can be extended. As mentioned in Sec. Ill, 
we considered the pairing interaction to the leading order in A///. The same approximation 
was made also in Refs. Since A/// can be as large as ~ 1/3, it is clearly 

desirable to go beyond this approximation. Secondly, the density dependence of A and m s 
can be included to get a more complete picture. (In the present work, we have varied A 
over a wide range and studied how it affects our results. It appears that allowing A to 
depend upon the density of the medium, may not change the results significantly. See, e.g., 
the two solid curves in Fig. 7(a), where we show the variation of the density as one goes 
from the center to the surface of the star. Each solid curve is for a fixed value of A, namely 
or 100 MeV. Obviously, curves for intermediate values of A (say, 10, 20, 90 MeV) 
would lie in between these two curves. Hence if A were allowed to vary with the density or 
equivalently with r, the resulting curve is not expected to be drastically different from the 
two solid curves shown here.) Thirdly, existence of additional phases (like 2SC, crystalline, 
etc.) can be looked into. Finally, we did not include perturbative corrections to the equation 
of state; these can be incorporated. Furthermore, a similar study where the bag parameter 
is dependent on density can be made for hybrid stars. 

We conclude by noting that the limitations of the present work, pointed out in the 
preceding paragraph, and the simplicity of the theoretical model employed, should be kept 
in mind while interpreting the numerical results presented here. 
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Figure Captions 



Fig. 1: Bag parameter B vs baryon number density p for the three models considered in 
this paper. 

Fig. 2: Energy density e, pressure P, bag parameter B, and \p dB/dp\ vs quark number 
chemical potential /i, for the gap parameter A = 100 MeV. Here and in the rest of the 
figures, the panels (a), (b), and (c) refer to the three models of the density-dependent bag 
parameter B(p), described in Sec. II. (For the sake of clarity, Bx5 rather than B is plotted.) 

Fig. 3: Pressure P vs energy density e, i.e., the equation of state of the quark matter. 
In each panel, the two solid curves correspond to a calculation with a density-dependent 
B, and the two dashed curves correspond to a calculation with a constant B. Curves are 
labelled by the values of A in MeV. 

Fig. 4: Mass M (in units of the solar mass M Q ) vs radius R of a strange star. Notation 
as in Fi gs. 2 and 3. The dot-dashed line is the M/M & ~ 0.15 R (in Km) constraint obtained 



Fig. 5: Mass M (in units of the solar mass M©) vs central energy density e c (in units of 
the energy density Eq of the normal nuclear matter) of a strange star. Notation as in Figs. 



Fig. 6: Radius R vs central energy density e c (in units of the energy density Eq of the 
normal nuclear matter) of a strange star. Notation as in Figs. 2 and 3. 

Fig. 7: Baryon number density p (in units of the baryon number density po of the normal 
nuclear matter) vs distance r, for the maximum-mass strange star. Notation as in Figs. 2 
and 3. 

Fig. 8: The maximum mass M of a strange star (in units of the solar mass M & ) vs the 
gap parameter A. Notation as in Figs. 2 and 3. 

Fig. 9: The maximum radius R of a strange star vs the gap parameter A. Notation as 
in Figs. 2 and 3. 
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